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Abstract 

The one-dimensional small-polaron model with open boundary conditions 
is considered in the framework of the quantum inverse scattering method. The 
spin model which is equivalent to the small-polaron model is the Heisenberg 
XX Z spin chain in an external magnetic field. For spin model, the solutions of 
the reflection equation (RE) and the dual RE are obtained. The eigenvalues 
and eigenvectors problems are solved by using the algebraic Bethe ansatz 
method. The case of the fermion model is also studied, the commuting of 
the transfer matrix can be proved, and the eigenvalues and the Bethe ansatz 
equations for fermion model are obtained. 
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1 Introduction 



Recently, much more attention has been paid to study the integrable models with 
open boundary conditions [1-10]. This was initiated by Sklyanin [1], who proposed a 
systematic approach to handle integrable models with reflecting boundary conditions 
in which the so-called RE plays a key role. 

Sklyanin assume that the R-matrix of the vertex models are P and T symmetric 
Ru = PRyiP = i?2i,-Ri2* 2 = R12, where t i} i = 1,2 means transposition in the 
i-th space, P is the permutation operator Px <g> y = y <g> x. Because only few 
models satisfy these properties, Mezincescu et al [3] extended Sklyanin's formalism 
for constructing integrable open chains to the case where the R matrix is only PT- 
invariant PR 12 P = R^f 2 - It is later realized that the unitarity and cross-unitarity 
relations of the R-matrix are enough to prove integrable open boundary conditions. 
In the framework of the quantum inverse scattering method (QISM) [11], a lot 
of models with open boundary conditions are solved by using the algebraic Bethe 
ansatz method [2-7]. And the boundary cross-unitarity is also proposed for open 
boundary conditions problems [8]. 

In this paper, we will study the one-dimensional small-polaron model which de- 
scribes the motion of an additional electron in a polar crystal [12]. Via the Jordan- 
Wigner transformation, the one-dimensional small-polaron model is equivalent to 
the Heisenberg XXZ spin chain with an external magnetic field parallel to the z 
direction [13]. We will consider the open boundary conditions for this Heisenberg 
XXZ spin chain. The R-matrix satisfy the unitarity and the cross-unitarity rela- 
tions, so we can formulate the RE and the dual RE. With the help of these relations 
we can prove the commuting of the transfer matrix with open boundary conditions. 
The commuting transfer matrix can give an infinite set of conserved quantities, it 
is general in this sense that the model under consideration is integrable. It is well 
known that the first non-trivial conserved quantities is the Hamiltonian. Solving 
the RE and the dual RE, we find the diagonal solutions K matrices of those REs. 
From the transfer matrix with open boundary conditions, we obtain the Hamiltonin 
which includes the boundary terms in the Heisenberg XXZ spin chain with an ex- 
ternal gnetic field parallel to the z direction. Those boundary terms are determined 
by the boundary K matrices. Using the algebraic Bethe ansatz method, we find 
the eigenvalues and eigenvectors of the transfer matrix. So the eigenvalues of the 
Hamiltonian can also be obtained. 

Using similar procedure, the one-dimensional small-polaron model with open 
boundary conditions is studied. The eigenvalues of the transfer matrix are obtained 
by using the algebraic Bethe ansatz method. 

The paper is organized as follows: In section 2, we will introduce the model and 
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formulate the RE and dual RE. We then prove the commuting of the transfer matrix 
with open boundary conditions in section 3. In section 4, we use the algebraic Bethe 
ansatz method to solve the eigenvalues and eigenvectors problems. Section 5 will 
be devoted to the fermion case, and the eigenvalues and the corresponding Bethe 
ansatz equations are obtained for the one-dimensional small-polaron model with 
open boundary conditions. Finally, we will give a summary and discussions. 



2 The model and the RE and dual RE 

As mentioned in the introduction, the one-dimensional small-polaron model de- 
scribes the motion of an additional electron in a polar crystal [12]. The Hamiltonian 
takes the form 

N N N 

H = Wj2 n j - J H( a ] a i+i + a )+i a i) + V H n o n 3+ii (!) 
j=l j=l j=l 

where Oj,cij are fermion creation and annihilation operators at lattice site j, and 
satisfy the anticommutation relations 

{di, aj} = {a\, a]} = 0, {dj, aU = (2) 



71; 



a\dj is the occupation number operator. W, J and V can be found in the paper 



3 ^3^3 



of Fedyanin et al [12]. 

Applying the the Jordan- Wigner transformation for a, ,aj and rij 

cij = exp ^7r^(7j + (7fj a~, 
a] = exp ^mt ^ <7, + of j o-/, 

n, = ^ (3) 

where af = \{a°- ± ia V j) and crj, a|, a| are Pauli operators at lattice site j. The 
Hamiltonian of the one-dimensional small-polaron model changes to the Heisenberg 
XX Z spin chain with an external magnetic field parallel to z direction. 

N / V \ w +v N 

H = -E [A^a- +1 + *t*+ +1 ) - a^ +1 ) + E a] + cost. (4) 

j=l ^ 4 / Z j=1 
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Here the model is considered in the periodic boundary conditions. 

From this Hamilotanian, the local monodromy matrix L operator is found to be 

[13] 



a',+y, a'.-b'. 



Lj(X) = 

here we have used the notations 
a' + : b' + : a'_ : b'_ : c = 

with 



+ 



CO, 



CO 



+ 



2 ®3 3 

a'_+b'_ a'_-b'_ 



2 -°3 



(5) 



£(A) sin(A + 2r/) : ^(A) sin(A) 

: r X (A) sin(A + 2?7) : f (A) sin(A) : sin(2r/), 



£(A) = sec (a) cos(A + a) sec(A). 

And J, V and W are given by 
V 

J : -— : (W + V) = 1 : cos(2r?) : 2 sin(2r?) tan(a). 



(6) 



(7) 



(8) 



It is well known that for periodic boundary condition problems, the Yang-Baxter 
relation [14,15] plays an essential role which takes the form 



i? 12 (A,/x)L 1 (A)L 2 (/x) = L 2 {ii)L 1 (\)R 12 (\n), 



(9) 



where the indices of L represent the auxiliary spaces, and the R-matrix can be found 
to be: 



#i 2 (A,aO = 



/o+(A,/i) \ 

6_(A,/x) c(A,/x) 

c(A,/x) b+(\,fi) 

\ a_(A,yu) / 



(10) 



the non-zero elements of the R-matrix are defined as 



a+(A, //) = £(A)r>) sin(A -fi-2rj), a_(A, //) = £ (A)£(/x) sin(A -fi-2 V ), 
b-(\, li) = C'WC 1 ^) sin(A - p), 6+(A, //) = £(A)£(/i) sin(A - ^), 

c(A,//) = sin(277), (11) 

here rj is the anisotropic parameter. In the framework of the QISM, The row-to-row 
monodromy matrix is defined as: 



T(A) = L JV (A)..-L 1 (A) 



;i2) 
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Here we know that the indices of the operator L mean the quantum spaces. It is 
known that using repeatedly the Yang-Baxter relation (9), we have 

R 12 (X,fi)T 1 (X)TM = T 2 (p J )T 1 (X)R 12 (X, f i). (13) 

For periodic boundary conditions, we define the transfer matrix as trace of mon- 
odromy matrix T in the auxiliary space. Here, to avoid confusion, we represent the 
auxiliary space by 0. 

tperi.W = ^qT^X) 

= troL JV o(A)---L l5 (A). (14) 
The Hamiltonian (4) can be obtained from this transfer matrix 

d 

H = -sin(2??)— \nt peri .(X)\x=o- (15) 

What we study in this paper is the reflecting open boundary conditions. For 
open boundary condition case, besides the Yang-Baxter relation, we also need the 
RE and the dual RE. Generally, we write the RE as the following form 

R 12 (X, f i)K 1 (X)R 21 ( f i,-X)KM = K 2 (p J )R 12 (X,- f x)K 1 (X)R 21 (- f x,-X). (16) 

The form of the dual RE depends on the unitarity and the cross-unit arity relations 
of the R-matrix. For the model considering in this paper, the unitarity and cross- 
unitarity relations of the R-matrix are written as 

i?i 2 (A, n)R 2 i(fJ,, A) = sin(2?? + A — //) sin(2r? - A + //) = p(X - //), (17) 
X%{\ p + 477)i&(/x, A) = - sin(A - fi) sin(A - /x - Arj) = p(\ - fi). (18) 

Here, for convenience, we have introduced the notations 

- Yl2<A ''' ) = {(^ ( -4^) Ml( '' )fll2(A ''' )Ml( ^' 

M ^-(l «A-4° r .(A))- < 19 ' 
It is also convenient to write here another form of the cross-unitarity relation: 

^ 2 (A, n)Y*{ (fi - A V , X)=p(X-n + A V ), 
Y 12 (X,fi) = ^ A + 4 ^ M 2 (A)i? 12 (A, j »)M 2 (A). (20) 



Thus, we take the dual RE as the following form: 

R 12 (-X, -fj,)Kf(\)Y 2 i(-fi - 4r7, mt^) 

= K+(p)X 12 (-X,p + 4r,)K 1 (X)R 21 (p,X). (21) 

Following the method of Sklyanin, we define the double-row-monodromy matrix as: 

T(A) = T{\)K{\)T-\-\). (22) 

Using the Yang-Baxter relation (13) and the RE, it can be proved that the double- 
row-monodromy matrix also satisfy the RE: 

R 12 (X, fi)T 1 (\)R 21 (n, -\)T 2 {n) = T 2 (fj,)R 12 (X, - fJ ,)T 1 (X)R 21 (-n, -A). (23) 

We define the transfer matrix for open boundary conditions as: 

t(X) = trK + {X)T{X). (24) 

3 Integrable open boundary conditions 

In this section, we will prove that the transfer matrix defined above for open bound- 
ary conditions constitute a commuting family. And we know that it is generally in 
this sense we mean the model is integrable. Now, let us prove the transfer matrix 
commute for different spectral parameters. 

t(X)t(fi) = tr 1 K 1 + (A)T 1 (A)tr 2 K 2 + (/i)T 2 (/i) 
= tr 12 K 1 + (A)* 1 K 2 + (^)T/HA)T 2 (^), 

here we take transpositon t\ in space 1. Now we insert the cross-unitrarity relation 
(18), we have 

tr 12 K?(X) h K+^)XU-X, ii + 477)i&(^ -A)T/ 1 (A)T 2 ( /U ), 

tr 12 {Kt{X)^Kt^)XU-\ P + 4r7)}* 1 {T 1 (A) J R 21 ( /i , -A)T 2 (/i)}, 

tr 12 {K+(p)X 12 (-X, ix + 4r^+(A)}{T 1 (A) J R 21 ( M , -A)T 2 (/i)}. 

In the above calculation, we take transposition in the first space. Insert the untarity 
relation of the R-matrix (17), and use the RE (23) and the dual RE (21), we have 

1 



1 


p(- 


-X- 






1 




p(- 


-A- 


p) 




1 




p(- 


-A- 





p(-X - n)p{X - n) 



tr 12 {K+(p)X 12 (-X, ii + A v )K+(X)R 21 (p, A)} 
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{R 12 (X,p)T 1 (X)R 21 (p,-X)T 2 (p)}. 

' -tr 12 {R 12 (-X, -p)K+(X)Y 21 (-p - A V , X)K+(p)} 



p(-X - p)p(X - p) 
{T 2 (p)R 12 (X, -p)T 1 (X)R 2 i(-p, -A)}. (25) 

Use again the unitarity relation (17) and the cross-unitarity relation (20), and con- 
sidering the properties of the function 

p(X-p) = p(p-X), 
p(-X-p) = p(X + p + Arj), (26) 



we have 

1 



-tr 12 {K+(X)Y 21 (-p - 477, A)K 2 + (^)}{T 2 (^)it> 12 (A, -^(A)}. 



p(-X- p) 

= p { _X-p) tri2{Y * l{ ~^ ~ 4?? ' A )^i + ( A )* 1 ^2 + (^)}{^(^)T/ 1 (A) J Ri 1 2 (A, -p)}. 
= tr 2 K+{p)T 2 {p)tr l Kt{X)T l {X) 

= t(p)t(X) (27) 

Thus we have proved that the transfer matrix t(A) forms a commuting family which 
gives an infinite set of conserved quantities. 

t(X)t(p) = t(p)t(X) (28) 

One solution of the RE and dual RE correspondes to one transfer matrix. We 
have seen that in the proof of the integrability, the RE and dual RE are independent 
of each other. 

By solving directly the RE and the dual RE, we find the following diagonal 
solutions to the RE and the dual RE: 

K (X) = ( n g(-A) sinftHA) J , (29) 
V U £(A) sin(V-A) / 

(I \ 

K ( A ) = n €W sm('4,+ ~\~2 V ) , (30) 
V U £(-A) sin(</>++A+2»j) / 

where if) and if) + are boundary parameters. We noticed that this diagonal solutions 
are also the general solutions of the RE (16) and the dual RE (21). 
Explicitly, the transfer matrix can be written as: 

f(A) = trM(X)L N - (X) ■ ■ • L l0 (A)Ko(A) Ao(A) • • • C N - (X). (31) 



Here we denote C(X) = L x (— A). Up to a whole factor, we take the inverse of the 
L operator as: 

/ + <^o* car \ 

A(A)= 2 + + 2 ] a » + «> l»- b » z • (32) 

where 

o^(A) : fe" (A) : a" (A) : b"_ (A) : c 
= r 1 (-A)sin(A + 27 ? ):£(-A)sin(A) 

: £(-A) sin(A + 2rj) : £ _1 (-A) sin(A) : sin(2r/), (33) 
We define the Hamiltonian for open boundary conditions H OB as follows: 

Hob = -^sin(277)-^lnt(A)| A=0 

Vtf Kin(<2n)K>(0) ^o^(0)^o(0)^o mi 
- - 2_, ^jJ+i - 2 sm (2^)^i(0) t r #+(o) ' ( ^ 

Here 

• i = L' jJ+1 (0)P jJ+1 = I ).,,.:£'{()). (35) 

And the relation 

tr- K^0)L' m (0)P m = tr- o K£(0)P N - o (0)jC' N - o (0) (36) 

has been used. 

So, the open Heisenberg XX Z spin chain with an external magnetic field parallel 
to z direction is as follows: 

N ~ 1 / y \ m/ + v N ~ 1 

Hob = - E [M'J+i + *7< + i) " + E -1 

i=l V 4 / Z j=1 

+i sin^)^^ + I (V + W - sin^)^^ + cost.. (37) 
2 sin -0 2 y sm -0+ J 

4 Algebraic Bethe ansatz method 

In this section, we will use the algebraic Bethe ansatz method to obtain the eigen- 
values and eigenvectors of the transfer matrix with open boundary conditions. As 
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mentioned above, the double-row-monodromy matrix is defined as T(A) = T(A) 
xi\f(A)T _1 (— A), and we denote it as follows: 

T <-» = (c(A) ^)=™A)7-(-A) 

/ -4(A) B(A) W 1 A / i(-A) fl(-A) \ 

" i 0(A) D(X) ) { ffigSl J ( C(-A) fl(-A) ) < 38 > 

We know that one can find easily the pseudovacuum state |0 >, and acting the 
elements of the monodromy matrix T(A) and T _1 (— A) on this pseudovacuum state, 
we have 

C(A)|0>=0, C(-A)|0>=0, 

5(A)|0>^0, B(-X)\0>^0, 

A(X)\0 >= a(X)\0 >, A(-X) =a(-A)|0 >, 

D(X)\0 >= S(X)\0 >, D(-A) = 5(-A)|0 >, (39) 



where 



a(A) = (e(A)sin(A + 2r ? )) iV , 

5(A) = (£(A)sin(A)f, 

«(-A) = (r 1 (-A)sin(A + 27 7 )) Ar , 

*(-A) = (r^-A)^))^. (40) 



From the Yang-Baxter relation (13), one can obtain 

T 2 - 1 ( f ,)R 12 (X, f i)T 1 (X)=T 1 (X)R 12 (X,n)T 2 -\ f ,). (41) 
Let pi = —A, use the result C(A)|0 >= 0, we have 

£(-A) sin(2,) D{x)D{ _ m> (42) 



f(A) sin(2A + 2??) 
For convenience, we introduce here a transformation 



Then for double-row-monodromy matrix T, acting its elements on the pseudovac- 
uum state, and considering the above transformation, we have 



C(A)|0>=0, £(A)|0>^0, 
.A(A)|0 >= a(A)a(-A)|0 >, 

f (-A) sin(2A) sin(^ + A + 2r?) 



V(\)\0 >-- 



f (A) sin(2A + 2rj) sin(^ - A) 



<J(A)5(-A)|0 > 



(44) 



With the help of the transformation (43), the transfer matrix for open boundary 
conditions can be written as: 



t(X) = A(\) + 



f(A) sin(^ + - A - 2rj) 



f (-A) sin(^+ + A + 2r]) 
sin(2A + Arj) sin(^+ + A) 



V{\) 



sin(2A + 2rj) sin(^+ + A + 2r]) 



£(A) sin(^ + - A - 2??) ~ , 



£(-A) sin(^+ + A + 277) ' 



(45) 



We know that the double-row-monodromy matrix T satisfy the RE, from the 
relation (23), we have the commutation relations: 



A(\)B(im) 



sin(A — // — 2rj) sin(A + fi) £ 2 (— A) 
sin(A - n) sin(A + Jt + 2r?) £ 2 (A) 
| sin(27y) sin(2//) ^(-A)e(-/i) 
sin(A — /i) sin(2/i + 277) £ 2 (/i) 
sin(2r7) £(-A) 



B(M)^(A) 
B(A)^) 



sin(A + ^ + 2?7) £(//) 



■B(A)2>(^) 



(46) 



P(A)B(|i) = 



B(A)B(/x) 



sin(A - /i + 2rj) sin(A + /i + 4r?) £ 2 (-A) 
sin(A - /j) sin(A + /i + 277) £ 2 (A) 
sin(2r?) sin(2A + 4r?) £ 2 (-A) 



sin(A - fi) sin(2A + 2r/) f(A)£(/i) 

sin(2r7) sin(2/j) sin(2A + 4r?) £ 2 (-A)£(-/i) 
+ sin(2A + 277) sin(2/j + 277) sin(A + Ji + 2r?) £(A)£ 2 (/j) 

*W(-a) 



ew(-p) 



B(X)A^), 
(47) 



Using the standard algebraic Bethe ansatz method, we assume the eigenvectors 
take the form 



•••,//„ >= • • •B(fi n )F(fji 1 ,- • -,//„)|0 >, 



(48) 
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where F is a non- vanishing function. Acting the transfer matrix on this eigenvector, 
we obtain the eigenvalues of the transfer matrix as: 



A(A,//i, • • •, (in) 

sin(2A + 477) sin(^ + + A) 
sin(2A + 2rf) sin(V>+ + A + 2r?) 

xa(A)a(-A) ft 1 2V) + "° ehX) 



i=i 



sin(A - /ij) sin(A + /i« + 2r/) £ 2 (A) j 
sin(^+ - A - 277) sin(2A) sin(^ + A + 2rj) ^ 
sin(^ + + A + 2r]) sin(2A + 2rj) sin(^ - A) 1 M J 
" [ sin(A + 2r,) sin(A + ^ + 4?]) £ 2 (-A) \ 
l\ \ sin(A - fj^) sin(A + ^ + 2»y) £ 2 (A) J ' 1 ' 

the parameters //j, i — 1, • ■ ■ ,n should satisfy the following Bethe ansatz equations: 

sin(/ij + 2i]) 2N sin(ip + — iij — 2rf) sin(^ + fij + 2ry) 
sin(/x,) 27V sin(-?/> + + fj,j) sin(^ — /ij) 

yr j sin(/Xj - //j + 277) sin(/ij + + 4t/) | 
iJJ^- 1 sin(^ - /ii - 2ry) sin(^ + J ' 

j = l,2,...,n. (50) 

The Bethe ansatz equations ensure that the unwanted terms vanish. On the other 
hand, the Bethe ansatz equations also ensure that the eigenvalues of the transfer 
matrix are entire functions. 

From the definition of the Hamiltonian (34), we thus can obtain the eigenvalues 
E of the Hamiltonian 

E = -Ncos(2r))-y — — — + (2n — N) sin(2?7) tan(a) 

- - . . — t . " — + cost. 51 

2sin(^+)sin(^ + + 2rj) v ; 

5 Integrable model for fermion case 

In this section, we will study the corresponding fermion model. The Lax represen- 
tation for the Hamiltonian (1) is found to be [13] 

L > W ~ \ -tea) a'_ - (a'_ + ib'_)n 3 ) ' (52) 
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Here a±,b'± and c have already be defined in the above sections. As for the spin 
model, we first write out the graded Yang-Baxter relation, 



i?f 2 (A, fi)L( {X)L F M = L F fr)L( (A)f?f 2 (A, p), 



(53) 



we should notice that the meaning of some notations here is different from that of 
spin models presented above. Here 



Lf (A) = L F (X) ® s 1, 



L F (X) 



L F (X). 



(54) 



Cg> s denotes the super tensor product 

[A ® s B] ijM = (-iy^ + ^A lk B 3l (55) 
with parity p(l) = 0, p(2) = 1. The R-matrix for the fermion model is defined as: 

\ 



( a+(X,n) 

—ib-(\,n) c(A,/x) 

c(A,/x) ib + (\, n) 

-a_(A,/x) 7 



(56) 



We can prove that this R-matrix satisfy the following unitarity and cross-unitarity 
relations 

i?f 2 (A, / u)i? 2 v 1 (/i, A) = sin(2r] + A — fi) sin(2r] — A + /i) = p(A — /i), (57) 
Xf^* 1 (A, /i + 477 — Tr)/?^ 5 * 1 (/!, A) = sin(A — /i) sin(A — /i — 477) = — p(A — /i), 



7? 



F s *i , 
12 



(A,/i)y 2 f* 1 ( / u-4r ? + 7r,A) = -p(A-// + 4?7), 



(58) 



where we have used the notations: 
x£(A,/x) 

M F (X) = 



f (/x - 4ry + tt) 



Mf(^i?f 2 (A,//)Mf(^, 



g(A+ ^ } 1t) m f {x)r f 2 {x^)m f {x). 

e(A-4r ? + 7r)r 1 (A) J' 



(59) 



Here, 



R 



21 



P F i?f 2 P F 



(60) 
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P F is super permutation operator which is defined as: 



PS,M = (-lf t)pU) Su6 jk (61) 
st means super transposition defined as 

We can also write the RE and dual RE for fermion models, here we take the solutions 
of RE and dual RE K F and K +F to be diagonal. It is trick to deal with non-diagonal 
K-matrix, because there are grassmann odd numbers in the non-digonal positions. 
One can write the graded RE and dual graded RE for fermion model as: 

R^(X,n)K[(X)R^(ji,-X)K[^) = K 2 (^) F J Rf 2 (A,-^)Kf(AX(-/i,-A). (63) 

i?f 2 (-A, -fi)Kf(\)Y 2 F (-v - 477 + 7T, X)Kf(n) 
= K+ F (»)Xf 2 (-\, fi + A V - tt)K[ (AX(/x, A). (64) 

By solving the graded RE and the dual graded RE, we find respectively the solutions 
of the graded RE and dual graded RE as: 

/ g(A) sin(V>-A) n \ 

K F (X) = I «-a)«jW+a) u j j (65) 

/ sin(V>++A+2r ? ) n \ 

K + (A) = f « A ) sin(V-+-A-2r,) U ^ j . ( 66 ) 

One can find K F is similar as the spin model, but there appear a — in K +F compared 
with spin model. The transfer matrix for fermion models is defined as: 

t F (\) = strK +F (X)T F (X) (67) 

Here, str means super trace defined as strA = (— l) p ^ A^. And 

T F {\) = L F (X) ■ ■ ■ L F (X)K F (X)£ F (X) ■ ■ ■ C F (X) (68) 

satisfy the graded RE. We have used the notation C F (X) = L F 1 (— A), which is 
defined as: 



C ' [ X) - I -ca) -o!'_ + {a'L - ib'L)n j I" 
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The definition for a±, b± and c can be found in the above sections. The Hamiltonian 
TioB for small-polaron model with open boundary conditions is as follows: 

N-l N-1 N-1 

Hob = W n j - J ( a i a i+! + a i+i%) + V n j n j+i, 
j=i j=i j=i 

(V . .cost/A (V . ,cosip + \ 

+ — + sin(277)- — - \n 1 + — + W - sm(2r))- — — n N + cost. 
\2 sin ip J \2 sm ip + J 

(70) 

The commuting of the transfer matrix for fermion model with open boundary 
conditions can be proved similarly as that for the spin model, here we should notice 
that in the proof of the integrability, besides the super transposition st, we also need 
a inverse of the super transposition st with A stst = A. 

We denote 

T F (X) - ( AF[X) EF[X) \ (71) 
y C F (X) V F (X) ) ' {71) 

thus we can write the transfer matrix explicitly as 

fF(x) C(-A) sin(^+ + A + 2rj) F F 
1 (A) " £(A) sin(^-A-2 ?? )' 4 {X)+V (A) 



e(-A)sin(^ + + A + 2r ? ) ~ F 



£(A) sin(^+ - A - 2rf) 



■A (A) 



sin(2A + 4r / )sin(^+-A) gF 
sin(2A + 2rf) sin(^+ - A - 2rj) 1 1 ' 

In above, as for the case of spin model, we have introduced the following transfor- 
mation for convenience, 

Following the standard algebraic Bethe ansatz method, we define the pseudovac- 
uum state for fermion model as |0 > F with a,j\0 > F = 0. One can find 

^(A)|0 > F = j(A) sin(2A) s y-A-2,) 

v 71 C( _ A) sin(2A + 2rf) sin(-0 + A) v 71 

D F (A)|0> F = 5 F (A)|0> F , 

£ F (A)|0> F =0, C F (A)|0> F ^0. (74) 
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Where we have 

a F (X) = (-^)^in^(A),^(A) = ("^) "^(A + 2rj). (75) 

Because the double- row-mono dromy matrix for fermion model satisfy the graded 
RE, we can obtain the commutation relations which are necessary for the algebraic 
Bethe ansatz method. 

a f (\)c"m = 3i " (A = " + f > 5 ;; 1(A + p> S} x I cma f w 

v ' K ' sm(A - fi) sm(A + fi + 2t/) £ 2 (-A) v ' v ' 

sin(27y) sin(2A + 4r/) e(A)£(/i) F ? F( , 
sin(A-/i)sin(2A + 27/) £ 2 (-A) { } {N 

sin(2r/) sin(2/i) sin(2A + 4r/) 
sin(2A + 2rj) sin(2/j + 2r}) sin(A + // + 2r}) 

x ^}fi^_ C F (x)v F (u) 

x e(-m-») c {x)v ^' 

v f {x)c f { ) = sin(A - /, - 2t/) sin(A + fl ) J 2 ^ F 
y ' ^' sin(A-/V)sm(A + /i + 277)£ 2 (-A) vry v 1 

sin(2r7) sin(2/x) £ 2 (/i) C F (\)V F (n) 



sin(A - //) sin(2/j + 27]) £(-A)£(-/j) 
sin(2?7) ^(/i) j, ~ F 

'sin(X + ^ + 2r ] )a-X) ^ { ^ 



(76) 



Those commutation relations are only slightly different from the relations for spin 
case. Assume the eigenvectors take the form C F (/ii) • • •C F (// ra )-F(/j 1 , • • • , /i n )|0 > F , 
where F(fj,i, ■ ■ • fj, n ) are non- vanishing functions. We can find the eigenvalues for the 
fermion transfer matrix with open boundary conditions are: 

A F (A,/^i, •••,/!„) 

sin(^ + + A + 2r]) sin(2A) sin(^ - A - 27/) F 
~ sin(^+ - A - 2t/) sin(2A + 2r/) sin(^ + A) " 

sin(A - Atj + 27/) sin(A + ^ + 4r/) £ 2 (A) 1 



x 



n 

i=i 



sin(A — Hi) sin(A + /Xj + 2rj) £ 2 (— A) 



sin(2A + Arj) sin(V;+ - A) ^ F 
sin(2A + 2rj) sin(^+ - A - 2rj) 
A f sin(A - ^ - 2r/) sin(A + ^) £ 2 (A) ] 
Hlsin(A-^)sin(A + /i l + 2r / )e 2 (-A)J • 1 J 
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The parameters /xi, • • • , fi n are restricted by the Bethe ansatz equations: 
sin 2Ar ( y u j + 2ri) sin(^ + + fij + 2rf) sin(^ - ^ - 2rj) 

yr j sin(/Xj -JM + 2rj) sin(^ +JM + 4rj) \ 
i=l ^ \ sin(^ -m- 2r]) sin(// j + J ' 

j = l,2,...,n. (78) 
The eigenvalue of the Hamiltonian for fermion model is: 

£ F = -ATcos(2?7) -V . + (2n — N) sin(2ry) tan(a) 

1 " ^ sin(^) sin(^ + 2??) 1 

sin 2 (2/7 ) 

+ cost. (79) 



2 sin(-0+) sin(-0 + — 2rj) 



6 Summary and discussions 

We formulated in this paper the RE and the dual RE for small-polaron model whose 
spin chain equivalent is the Heisenberg XX Z chain with an external magnetic field 
parallel to z direction. We found solutions to the RE and dual RE for both spin and 
fermion models. Using the algebraic Bethe ansatz method, we find the eigenvalues 
of the transfer matrix with open boundary conditions. 

It is interesting to use some results of this paper to calculate some physical 
quantities, such as the surface free ennergy and finite-size corrections, by using the 
thermodynamic Bethe ansatz method. The Lax pair method for open boundary 
conditions for this model is also an interesting problem. 
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